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CO 

Binarity has been hypothesised to play an important, if not ubiquitous, role in the formation of planetary nebulae (PNe). Yet there 
remains a severe paucity of known binary central stars required to test the binary hypothesis and to place strong constraints on the 
^ . physics of the common-envelope (CE) phase of binary stellar evolution. Large photometric surveys offer an unrivalled opportunity to 

efficiently discover many binary central stars. We have combined photometry from the OGLE microlensing survey with the largest 
sample of PNe towards the Galactic Bulge to systematically search for new binaries. A total of 21 periodic binaries were found 
thereby more than doubling the known sample. The orbital period distribution was found to be best described by CE population 
synthesis models when no correlation between primary and secondary masses is assumed for the initial mass ratio distribution. A 
comparison with post-CE white dwarf binaries indicates both distributions are representative of the true post-CE period distribution 
with most binaries exhibiting periods less than one day. An estimated close binary fraction of 12-21% is derived and is the first 
robust and independent validation of the prevailing 10-15% fraction estimated by Bond (2000). This suggests that binarity is not a 
precondition for the formation of PNe and that close binaries do not play a dominant role in the shaping of nebular morphologies. 
Systematic effects and biases of the survey are discussed with implications for future photometric surveys. 
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\ 1, Introduction observed period distributions of post-CE binaries that are very 

sensitive to the assumed values of acE and the initial mass ratio 

The common-envelope (CE) phase of binary stellar evolution de- distribution (IMRD) of binaries before the CE phase. The IMRD 

K> scribes a close binary immersed within a large envelope of ma- itself is a highly contentious field of study subject to many se- 

rN , terial lost from the primary component (Iben & Livio 1993). It i ection effects (Trimble 1990; Goldberg, Mazeh & Latham 2003; 

03 is established after mass transferred to the secondary overflows Halbwachs et al. 2003). Comparison of the period distributions 

the Roche lobe surrounding the secondary and engulfs the binary of known post-CE binaries therefore offers a valuable avenue to 

system. Energy and angular momentum are transferred from the constrain acE and the IMRD. 
binary system to the CE and the components are forced closer 

together. Evidence for the CE phase is seen in a diverse range The short ~ 10 y ear lifetime of P^netary nebulae (PNe) 

of short orbital period phenomena that include cataclysmic vari- means that close binar y CSPN are ldeal P robes °f CE evolu ~ 

ables, subdwarf B or white dwarf binaries, low-mass X-ray bina- tion because the orbital period distribution is observed directly 

ries, novae, degenerate binaries such as type la supernovae pro- a ^ ^ e end °f ^ fie phase where there has been no chance for 

genitors and the close binary central stars of planetary nebulae significant orbital evolution to take place as in older systems like 

(CSPN) cataclysmic variables. However, this potential has yet to be fully 

^ . . • i i • u- w ^ realised due to the previously very small sample size of binary 

Despite its crucial role in binary stellar evolution, there are ^gp^ 

critical aspects of the CE phase that are not well understood. 

This is exemplified by considerable uncertainty in the value of De Marco, Hillwig & Smith (2008, hereafter DM08) sum- 

the ocqe parameter that measures the conversion efficiency of or- marise the known close binary CSPN population of around a 

bital energy required to eject the envelope. Numerous population dozen objects. The rather inhomogeneous sample results mainly 

synthesis models have explored the role of acE albeit with vary- from the photometric monitoring of -100 CSPN over 30 years 

ing model assumptions and definitions of acE (de Kool 1990, after many different observing campaigns by Bond and collabo- 

1992; de Kool & Ritter 1993; Yungelson, Tutukov & Livio 1993; rators (Bond et al. 2000 and ref. therein, hereafter BOO). The sur- 

Han et al. 1995; Politano & Weiler 2007). The models predict vey biases and characteristics are therefore not well understood 
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(DM08) and these are compounded by a somewhat necessary 
bias towards low surface brightness PNe to reduce the effects of 
nebular contamination. The very small sample size and reserva- 
tions concerning the survey biases have hindered comparisons 
with predicted period distributions required to constrain acE or 
the IMRD. 

BOO estimate a close binary fraction of 10-15%, however 
DM08 expressed concerns about the accuracy of this very impor- 
tant quantity because of the survey biases. The close binary frac- 
tion is incredibly important for two interdependent reasons: (i) 
To investigate theoretical claims that PNe form largely, if not en- 
tirely, as a result of binary evolution (Paczynski 1985; Moe & De 
Marco 2006), and (ii) To determine the degree of influence bi- 
narity has in shaping the varied and complex nebular morpholo- 
gies of PNe (e.g. Balick & Frank 2002; Soker 1997). The latter 
topic has been the subject of considerable debate for the last 30 
years. Binarity offers a simpler or even preferred mechanism to 
explain non- spherical nebular morphologies such as bipolar PNe 
(e.g. Soker 1998), but without a firm estimate of the binary frac- 
tion other mechanisms cannot be ruled out (for a review of these 
other mechanisms see Balick & Frank 2002). 

New surveys for close binary CSPN are therefore required 
to make substantial progress in these areas by properly charac- 
terising the population and obtaining an independent estimate 
of the binary fraction. Ideally a photometric and radial velocity 
(RV) survey of all visible central stars of the -200 local volume 
PNe identified within 1-2 kpc (Frew 2008) would be conducted. 
However, this is difficult work for such a sparsely distributed 
sample and survey biases may be difficult to control (e.g. sam- 
pling and weather conditions). The most efficient and produc- 
tive, albeit magnitude limited, would be a photometric survey of 
a large population of PNe either towards the Galactic Bulge or 
Magellanic Clouds. Fortunately, the multi-epoch photometry for 
such an endeavour already exists thanks to microlensing surveys 
dedicated to these regions. 

After comprehensive exploration of the 
SuperCOSMOS Ua Survey (SHS, Parker et al. 2005) the 
Macquarie/AAO/Strasbourg Ha PNe catalogues (Parker et 
al. 2006; Miszalski et al. 2008a) have more than doubled the 
number of PNe towards the Galactic Bulge. We have combined 
these new discoveries with already catalogued PNe to search 
for new binary CSPN towards the Bulge using photometry with 
excellent spatial and temporal coverage from OGLE-III (OGLE; 
Sec. 12.11) . This survey dwarfs all previous efforts by examining 
nearly 300 PNe in an efficient and relatively uniform manner. 
Miszalski et al. 2008b presented some initial findings of the 
survey which we describe in full here. 

Section [2] introduces the observational data used during the 
search method described in Sec. [3] We present the new binaries 
and compare the period distribution with population synthesis 
model predictions in Sec. |U An estimate of the binary fraction 
is also made in Sec. [Hand we conclude in Sec.[5j 

2. Observations 

2.1. OGLE photometry 

The Optical Gravitational Lensing Experiment (OGLE) consti- 
tutes an extensive /-band photometric database constructed over 
many years of monitoring the Galactic Bulge and Magellanic 
Clouds. OGLE-I, being the first phase of the experiment, moni- 
tored a few select fields towards the Bulge (Udalski et al. 1992). 
The second phase, OGLE-II, brought much improved tempo- 
ral and spatial coverage of the Galactic Bulge (Udalski et al. 



2002b;Wozniak et al. 2002), but it was not until the third and 
most recent phase, OGLE-III (Udalski et al. 2002a, 2008) that 
OGLE photometry became suitable to uniformly assess CSPN 
binarity for a large population of PNe. OGLE-III is very well 
suited to the task with both large spatial coverage closely match- 
ing the spatial distribution of PNe (Fig. [I]), and well sampled 
fields sensitive to detecting the short periods of binary CSPN. 
Indeed, OGLE-III is effectively a 'one- stop shop' that super- 
sedes previous piecemeal and biased photometric CSPN mon- 
itoring campaigns from which any binary fraction estimates 
are rather uncertain. We therefore focus almost exclusively on 
OGLE-III in this work, but we have made some use of OGLE-II 
data in a reappraisal of Miszalski et al. (2009) to reveal one new 
binary (see later). 

The OGLE-III survey uses the 1 .3-m Warsaw telescope at the 
Las Campanas Observatory, Chile. It is equipped with a CCD 
mosaic camera with eight 2K x 4K CCDs giving a 35 x 35.5 
arcmin 2 field of view with 0.26" pixels sampling a median see- 
ing of 1 .2" (at worst 1 .8"). There are 267 survey fields tiled such 
that regions of highest stellar density and lowest reddening are 
favoured, but 80 of these have not yet been observed (Fig. [T]). 
Fields with higher stellar density are sampled more frequently. 

A non-standard /-band filter is used with a bandpass that ex- 
tends past the standard -9000 A drop-off (Udalski et al. 2002b) 
and reaches a limiting magnitude of / ~ 20. Nebular contami- 
nation can occur for PNe with strong [SIII] /19069, /19532 emis- 
sion lines that can be bright even in very reddened Bulge PNe 
(e.g. Jacoby & Van de Steene 2004). Other prominent lines such 
as [Arlll] 7750, [CI IV] 8045, Hell 8237 and members of the 
Paschen series lie favourably within the filter band-pass and 
could make a significant contribution if the reddening is not too 
high. A strong nebular continuum may also influence the /-band 
image in the brightest PNe. 

The vast amounts of OGLE-III data are processed by a data 
reduction pipeline based on the difference image analysis (DIA) 
method (Alard & Lupton 1998; Wozniak 2000; Udalski 2003). 
DIA essentially measures intensity differences between an ob- 
served frame and a high quality reference image to produce bet- 
ter quality time series photometry with lower scatter than other 
methods (Wozniak 2000). Systematic effects are generally very 
small but can become problematic for PNe where nebular con- 
tamination is high (Section l34l) . 

2.2. Optical Spectroscopy 

To complement OGLE-III photometry we have procured optical 
spectroscopy for a large number of PNe towards the Bulge with 
the 2dF/AAOmega and FLAMES multi-object spectroscopy fa- 
cilities. Our extensive spectroscopic data goes deeper than com- 
parable work conducted on 2-m class telescopes (Cuisinier et 
al. 2000; Gorny et al. 2004; Exter, Barlow & Walton 2004) 
enabling detection of many faint diagnostic features which we 
use to eliminate contaminants from our PNe sample (Sec. 13.51) . 
FigureQ]depicts the spatial coverage of AAOmega (large circles) 
and FLAMES (small circles) fields that overlap a large fraction 
of the better sampled OGLE-III fields. Additional spectroscopy 
is sourced from the MASH project (Parker et al. 2006; Miszalski 
et al. 2008a). The spectroscopic data are described below. 

The 2dF multi-object spectroscopy facility on the 4-m 
Anglo- Australian Telescope enables 400 objects to be observed 
simultaneously over a 2-degree field of view using optical fibres 
(Lewis et al. 2002). The fibres are positioned in the focal plane 
robotically and feed the AAOmega optical spectrograph (Sharp 
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Galactic Longitude 

Fig. 1. Galactic distribution of 35.5 x 35 arcmin 2 OGLE-III fields (squares) and PNe (points). Fields are coloured based on the 
number of total observations N, with warmer colours indicative of greater N (see key for N range). AAOmega fields (large circles) 
and FLAMES fields (small circles) indicate regions with relatively optical deep spectroscopy. 



et al. 2006). A total of 12 AAOmega fields were observed on 
our behalf in service mode on 26-27 March 2007 and 16 March, 
29-30 May and 8 August 2008. We used the 580V and 385R 
volume-phase holographic gratings to give central resolutions 
of 3.5 A and 5.3 A (FWHM) in the blue and red spectrograph 
arms, respectively. The standard dichroic was used to provide a 
wavelength coverage of 3700-8850 A, however the August field 
(€ = -2.5, b = -4.3) was observed with the redder dichroic 
which covers 4800-9740 A. Exposure times ranged from 2-3 x 
1200 s to 1-2 x 1800 s and shorter exposures were taken to en- 
able measurement of bright lines. Data were reduced using the 
2dfdr pipeline provided by the A AO, however the splicing of 
blue and red arms was performed separately. 

AAOmega field centres were largely selected to maximise 
target density, but some centres were modified to include bi- 
naries discovered during the course of this project. For 10 of 
the 12 fields we carefully inspected full resolution SHS and 
SuperCOSMOS Sky Survey (SSS, Hambly et al. 2001) data us- 
ing the techniques described by Miszalski et al. (2008a) to select 
fibre positions for known, MASH and MASH-II PNe. Many un- 
catalogued objects exhibiting Ha emission were assigned spare 
fibres and this eventuated in ~40 new very low surface bright- 
ness MASH-II PNe. These discoveries contributed to our master 
catalogue (Sec. 13. II) and some were covered by OGLE-III. This 
app roach also ensured many PN-mimics were observed (Sec. 
13. 51) . In most fields a few hundred sky targets were included to 
assist subtraction of the variable sky background. Prior to ob- 
servation fibres were assigned to selected positions using the 
simulated annealing field configuration algorithm (Miszalski et 
al. 2006). The relatively low target densities per field (-20-70 
PNe) enabled -200 sky targets and uncatalogued Ha emitters to 
be allocated, although the latter were given lower priority than 
catalogued PNe. 

Additional spectroscopic observations were conducted by 
AA and BM using the VLT FLAMES multi-object spectroscopy 
facility (Pasquini et al. 2002) during the ESO visitor mode pro- 
gram 0.79.D-0764(A) on 9-12 June 2007. The 25 arcmin diam- 
eter field of view and inclement weather limited the observed 



sample size. We used the deployable miniature integral-field 
units (mini-IFUs) that feed 20 optical fibres to sample a 2 x 3 
arcsec 2 field of view for each PN within a field. In this work we 
only make use of one field (1 = 4,/? = -3) that happened to con- 
tain three PN-mimics (Sec. 13.51) . Here we only make use of the 
30 and 60 minute exposures made with the GIRAFFE LR02 and 
LR03 filters, respectively, that gave a wavelength coverage of 
3948-5074 A. Data were reduced using the IRAF task dofibers 
and relative flux calibration was achieved using the spectropho- 
tometric standard star Feige 66. 

In the absence of AAOmega or FLAMES observations we 
draw upon the extensive spectroscopic resource of the MASH 
project (Parker et al. 2006; Miszalski et al. 2008a). The pre- 
dominantly confirmatory MASH spectra are not as deep, but are 
nonetheless useful for bright emission line ratios and in some 
cases detecting PN-mimics. Of most relevance to this work is the 
6-15 May 2008 MASH observing run on the Australian National 
University 2.3-m telescope. While the focus was to measure ra- 
dial velocities of a few hundred PNe towards the Bulge via short 
exposures (Miszalski et al. in prep), a number of the binaries dis- 
covered during the course of this project were given special at- 
tention. We used the 1200B and 1200R gratings with a slit width 
of 2" to give wavelength coverage windows of 4030-5050 A 
and 6245-7250 A at a resolution of -1.6 A (FWHM). Data re- 
duction was performed using IRAF with the assistance of the 
PNDR package!] 

3. Search Method 

Our search method involves a number of stages designed to find 
variable CSPN while incrementally refining the sample of PNe 
considered in later binary fraction calculations. After construct- 
ing a catalogue of PNe covered by OGLE-III (Section[3J]), iden- 
tifications of genuine or candidate CSPN are attempted using 
available near-infrared and optical images (Section EOl) . We then 
discuss the influence of nebular contamination (Section [33l) and 

1 http ://w w w. aao . gov. au/local/w w w/brent/pndr 
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general systematic effects (Section IX4l) on time series photome- 
try extracted for all identifications. Combining photometry with 
available spec trosc opy we proceed to remove non-PN contami- 
nants (Section [33]), before searching for both periodic (Section 
13.61) and aperiodic (Section IT71) variables indicative of binarity. 

3.1. Planetary Nebulae Catalogue Construction 

Full exploitation of OGLE-III for this work requires a large sam- 
ple of PNe towards the Galactic Bulge that must be compiled 
from various extant catalogues. We have therefore carefully con- 
structed a comprehensive PNe sample for \€\ < 13 and \b\ < 7 to 
cover all OGLE-III fields with coordinates of each PN verified 
using the SHS. The Macquarie/AAO/Strasbourg Ha PNe cata- 
logues MASH-I (Parker et al. 2006) and MASH-II (Miszalski 
et al. 2008a) form the foundation of the catalogue with -380 
and -100 PNe located in the region, respectively. An additional 
-450 PNe were then added from a large number of literature 
sources (e.g. Acker et al. 1992, 1996; Kohoutek 1994, 2002; 
Zanin et al. 1997; Beaulieu, Dopita & Freeman 1999; Cappellaro 
et al. 2001; Jacoby & Van de Steene 2004; Boumis et al. 2003, 
2006; Gorny 2006). As a further check we made use of updated 
coordinates provided by Kerber et al. (2003) and SIMBAD to 
ensure no objects were missed. 

Figure [T] depicts the most complete catalogue to date of 
Bulge PNe overlaid on OGLE-III fields. The OGLE field findeiQ 
reported that 297 out of the -930 PNe were covered by OGLE- 
III fields with at least one observation. These numbers do not 
include non-PN contaminants (Sec. 13.5b . The catalogue ex- 
hibits considerable variety from extremely low surface bright- 
ness MASH PNe to very high surface brightness non-MASH 
PNe such as NGC 6565. 

3.2. Central Star Identification 

Table IA.1I presents our sample of 297 PNe for which we have 
carefully examined all available images to attempt CSPN identi- 
fication. Central to our search was the inspection of 2x2 arcmin 2 
OGLE-III /-band images. The location and extent of each PN 
was established using an optical colour-composite image derived 
from the SHS and SSS with red, green and blue channels con- 
structed from Ha, Short Red and Bj images, respectively (e.g. 
Miszalski et al. 2008a). Visual comparisons between the opti- 
cal and /-band images were made using a near infra-red (NIR) 
colour-composite image as a reference point with red, green and 
blue channels constructed from the 2MASS K s , H and / at- 
las images, respectively (Skrutskie et al. 2006). OGLE-III V- 
band images and magnitudes were unavailable given the provi- 
sional calibration of Bulge data (Udalski et al. 2008), but for four 
PNe we have inspected their V-band images (effectively shallow 
[OIII] images). 

Each entry in Tab. IA.ll is remarked according to the type and 
quality of each identification (parenthesised and quoted labels 
discussed below). A 'true' (T) classification is based on either 
the presence of a single star at the geometric centre of the nebula 
or blue colours. Without the V-band data we could not search 
all PNe for blue CSPN with intrinsic colour (V - /) = -0.40 
(Ciardullo et al. 1999). Instead we have constructed SSS Bj-R F 
images that are sensitive in some cases to faint blue (B) CSPN 
(Parker et al. 2006). The 2MASS colour-composite is also useful 
in which bright CSPN of typically non-MASH PNe appear vio- 
let. Some blue CSPN were found differently, e.g. ESO Imaging 

2 http://ogle.astrouw.edu.pl/radec2field.html 



Survey data (Nonino et al. 1999) were used for PHR 1753-3443, 
and OGLE-III V-band images were used for H 2-29 and K 6-34. 
There are 92 PNe with 'true' CSPN identifications. 

In more ambiguous cases where a number of CSPN candi- 
dates are present we can only assign 'likely' (L) or 'possible' (P) 
based on available data. For a PN marked as 'likely' a number of 
CSPN candidates can be found near the geometric centre of the 
PN. If 'possible', then these candidates may not be near the geo- 
metric centre or the centre may be poorly defined. In both cases 
time series photometry is extracted and examined for all candi- 
dates bounded by the Ha emission of each PN. In practice dis- 
tinction between 'likely' and 'possible' objects is rather difficult 
as it is complicated by extinction and crowding. Unfortunately, 
the poor resolution of the often- saturated SSS data means few 
meaningful colour selections can be made resulting in the num- 
ber of 'true' CSPN identifications being a lower limit until V- 
band data become available. However, in some cases comparison 
of on- and off-band images from Ruffle et al. (2004) helped to 
better establish the geometric centre enough to select one CSPN 
as 'true'. There are 102 and 48 PNe with 'likely' and 'possible' 
CSPN identifications. 

A number of different scenarios occur when no CSPN iden- 
tification is possible. When the /-band image is available and 
absolutely no CSPN candidates were detected these are marked 
as 'NONE' (2 PNe). A large group of 34 'NULL' PNe include 
very bright non-MASH PNe avoided by the OGLE-III pipeline, 
PNe that fell just outside the requested /-band image or near a 
bad CCD column, and PNe for which no time series photometry 
could be extracted by the pipeline (Sec. 13.41) . Nebular detections 
are indicated by 'NEB' for 10 PNe only when there is a nebu- 
lar rim or other faint detection not coincident with a CSPN (e.g. 
NGC 6565, M 1-29). No label is given for all nebular detections 
but those indicated by 'NEB' are in the minority (see Section 
13.31) . Suspiciously spurious lightcurves caused by systematic ef- 
fects are indicated by 'S' for 14 PNe (see Sec. 13.41) . Periodic 
variables are indicated by 'BIN' (binary) and aperiodic variables 
by 'BIN?' (suspected binary). 

3.3. Nebular contamination 

Nebular detections occur for ~20% of the PNe in our sample to 
varying degrees due to the strong nebular [SIII] emission lines 
passed by the broad OGLE-III /-band filter (see Sec. 12. lb . Figure 
[2 shows a sample of nebular detection levels in our sample. 
Notably absent from Fig. [2] are MASH and MASH-II PNe of 
which almost none have nebular detections because they are in- 
trinsically much fainter than non-MASH PNe. 

In many nebular detections the nebulae appear faint or dif- 
fuse (e.g. K 5-20) and prove to be no obstacle to the stellar- 
optimised OGLE-III pipeline. Problems can occur however 
when the nebula becomes bright and large enough to appear non- 
stellar. Fully resolved nebulae such as M 1-29 and NGC 6565 
cannot be treated in any reasonable fashion (remarked 'NEB' in 
Tab. IA. lb . These nebulae may also confuse the pipeline even if 
a CSPN is clearly visible (e.g. Hb 4). Between these two ex- 
tremes are objects such as M 2-8 and M 3-14 where a bright and 
just resolved nebular core is present. The non- stellar appearance 
of these objects means extraction of time series photometry by 
the pipeline is not guaranteed. If photometry is obtained then we 
assume the nebula is obscuring an underlying CSPN and that 
we are sensitive to any variations present. This is a reasonable 
assumption given the brightness of the objects and indeed we 
have found one binary belonging to this description (H 1-33, see 
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Fig. 2. A selection of OGLE-III /-band images of PNe with neb- 
ular detections. All images are 30 x 30 arcsec 2 with North to top 
and East to left. 



later). These objects may however be more prone to systematic 
effects as discussed in the next section. 



3.4. Systematic effects 

The time series photometry of -14 PNe appear to be dominated 
by systematic effects producing spurious periodicities of one 
year or non-repeating 'outbursts' (Fig. [3]). While systematic ef- 
fects in OGLE-III are generally very small permitting transiting 
planet discovery, they can become prominent and difficult to ex- 
plain in extreme non- stellar objects such as PNe that show strong 
nebular contamination. In the most extreme PNe, i.e. many of 
those marked as 'NULL' , no time series photometry could be se- 
cured by the stellar-optimised pipeline (Wozniak 2000). In prac- 
tice only a few PNe may be rejected based on strong nebular 
contamination alone. The majority of 'NULL' PNe can rather be 
attributed to a very bright object being avoided by the pipeline 
or a bad CCD column across the PN in a small number of cases. 
Those that do pass the pipeline with significant nebular contami- 
nation may be more prone to systematic effects. Wozniak (2002) 
estimated 10% of OGLE-II Bulge variables are spurious hav- 
ing been caused by problems undetected at the pipeline level 
(Wozniak 2000). Some documented sources of systematic errors 
include colour effects (Udalski et al. 2002b) and 'mirrored' vari- 
ability from bright variables up to 100 pixels away (Mizerski & 
Bejger2002). 

An exact cause of the suspect variability in Fig. [3] is unclear 
based on the available data. The most likely cuplrit is differen- 
tial atmospheric refraction (DAR) exacerbated by strong nebular 
contamination. DAR is the shift in centroid position of a star de- 
pendent on the colour of the star and the airmass of the observa- 
tion (Alcock et al. 1999a, 1999b). When the nebula point- spread 
function is considerably larger than regular stars we suspect a 
small shift in the real centroid, compared to the fixed reference 
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Fig. 3. Examples of suspected spurious time series photometry. 
H 1-58 and M 2-23 show distinctly regular one year periods. Bl 
O and PHR 1752-3330 (prefixed by OGLE-II data) show unusu- 
ally long sustained 'outbursts'. 



image, to cause quite large light variations. Indeed, the 'inverted- 
U' shape present in H 1-58 and M 2-23 appears related to air- 
mass changes as the position of the Bulge in the sky changes 
during each observing season. Similar features are mentioned 
by Alcock et al. (1999a) as 'slow seasonal rolls' with one year 
periods. 

Photometric outbursts may be related to transient mass-loss 
events as in Lo 4 (Bond 2008), but we discount this explanation 
for Bl O and PHR 1752-3330 given the 1-2 year non-repeating 
duration is considerably longer than Lo 4 which is active ~8% of 
the time (Bond 2008). This is supported by inspection of differ- 
ent time series obtained at opposing edges of the fully resolved 
nebula of Bl O (Ruffle et al. 2004) where we found the two 'out- 
burst' peaks to suspiciously halve in magnitude alternately on 
each edge. 

Another source of systematic error pertains to those light 
curves that show considerably large scatter cr despite relatively 
small error bars of individual data points. This may be attributed 
to non-photometric or cloudy nights rather than large intrinsic 
variability (Udalski et al. 2008). We cleaned our sample by first 
selecting objects with large cr as a function of average magnitude 
(amongst all identified CSPN), and then checking whether the 
variability persisted when compared to stars surrounding each 
object. This process showed the vast majority of these cases were 
not caused by intrinsic variability (see also Sec. 13.71) . 

At this stage it suffices to identify any time series photometry 
with suspicious one year periodicity or outbursts. These have 
been remarked with an 'S' or 'S?' in Tab. IA.ll We refrain from 
marking PHR 1752-3330 as such given the outburst occurs for 
OGLE-II data only. Further investigation into systematic effects 
may be possible once V-band images and magnitudes become 
available to investigate the role of colour in DAR. 

3.5. Removing PN Mimics 

A significant number of PN mimics are present in existing PN 
catalogues (for an excellent review see Frew 2008; see also 
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Parker et al. 2006, Miszalski et al. 2008a). Mimics are often 
mistaken for PN based solely on their Ha emission without any 
deep follow-up spectroscopy to confirm their nature. We can- 
not overemphasise the importance of removing mimics to both 
reduce the possibility of identifying non-PN binaries and to ob- 
tain an accurate estimate of the binary fraction. We were able to 
identify and remove mimics based on 2MASS colours and deep 
FLAMES and AAOmega spectroscopy of identified Ha emitters 
in many fields overlapping with highly sampled OGLE-III fields 
(Fig.Q]). Table Q] contains the mimics identified in OGLE-II and 
OGLE-III fields that are discussed in the following sections. 





3.5.1. Symbiotic Stars 

In a photometric search for binary CSPN the main source of 
contaminants are symbiotic stars. Symbiotic stars are binary 
systems typically composed of a late-type giant or Mira trans- 
ferring material to a hot white dwarf companion via a stellar 
wind (Mikolajewska 2003, 2007) and they exhibit large and 
varied photometric variability (Mikolajewska 2001; Gromadzki, 
Mikolajewska & Borawska et al. 2007). Many are associated 
with extended nebulae (Corradi et al. 1999, 1995) that are easily 
mistaken for PN but are not the ejecta of the hot component as 
in PNe (Corradi 2003). Their PN-like nebulae can even exhibit 
peculiar variability (Doyle et al. 2000; Shaw et al. 2007). 

NIR colours can be very useful in identifying symbiotics 
(Schmeja & Kimeswenger 2001; Phillips 2007; Corradi et al. 
2008), but ultimately deep spectroscopy is required to identify 
typically faint features needed for definitive classification. These 
features include Raman- scattered OVI emission lines XX 6830, 
7082 unique to symbiotics (Schmid 1989), [OIII] /14363 useful 
in separating PN from symbiotics in the Rl = [OIII] /14363/Hy, 
R2 = [OIII] X5007/HJ3 plane (Gutierrez-Moreno, Moreno & 
Cortes 1995) and stellar features typical of the late-type compan- 
ion (Munari & Zwitter 2002; Belczynski et al. 2000). Indicative 
values of Rl ^ 0.6 and R2 ^ 10 are a reasonable signature of 
a symbiotic star but additional criteria should be satisfied for a 
definitive classification. 

During the catalogue construction phase (Sec 13.1b well- 
known symbiotics were removed using SIMBAD and litera- 
ture sources (Belczynski et al. 2000; Phillips 2007; Schmeja 
& Kimeswenger 2001). Application of the above criteria to the 
penultimate catalogue produced the cleaned final sample of 297 
PNe (Tab. I A. ID . We found the depth and large wavelength cover- 
age of our AAOmega data especially helpful in identifying sym- 
biotics. Almost all objects in Tab.[T]are newly confirmed (Sy) or 
strongly suspected (Sy?) symbiotic stars and many were previ- 
ously classified as PNe. Those presented here for the first time, 
with this work as the reference, were identified during the con- 
struction of the MASH (PHR objects) or MASH-II (MPA ob- 
jects) catalogues. M 2-29 has been classified as a PN harbouring 
a triple system (Hajduk et al. 2008), but we suspect a symbi- 
otic star after our spatially resolved FLAMES spectra confirmed 
high Rl and low R2 values near the central star (Fig. [4]) that were 
first reported by Torres-Peimbert et al. (1997). However, there 
may be another explanation for these line ratios and our tenta- 
tive classification will be refined in a future paper (Miszalski et 
al., in prep). Figure \5\ presents a sample of the extreme variabil- 
ity seen in these new symbiotics typified by irregular and long 
term variations. Outside OGLE fields we have reclassified more 
PNe as symbiotics that will be presented elsewhere along with 
more detailed analysis of the symbiotics in Tab. [TJ 
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Fig. 4. Spatial distribution of M 2-29 emission line ratios Rl 
(left) and R2 (right) suggestive of a symbiotic star. A star marks 
the approximate central star position that shifts in each exposure 
because of differential atmospheric diffraction. The FLAMES 
mini-IFU field of view is 2 x 3 arcsec 2 with North to top and 
East to left. 
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Fig. 5. A selection of OGLE-III time series photometry for con- 
firmed and strongly suspected symbiotic stars removed from our 
PNe sample. 



3.5.2. PHR 1805-2659: An M-dwarf flare star 

We reclassify the 'possible' MASH PN PHR 1805-2659 as an 
M-dwarf flare star (e.g. Doyle et al. 1990). Figure [6l shows the 
1800 s AAOmega spectrum taken on 29 May 2008 and OGLE- 
III time series photometry of the object. The M-dwarf spectrum 
exhibits strong Hydrogen B aimer and Ca H and K emission lines 
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Remarks: 

Sy(?): Symbiotic star (suspected); F*: Flare star; U: Unlikely PN; A: AAOmega observations 

F: [FeVII] /U5721, 6087 present; H: Broad Ha; M: Late-type features; N: Resolved Nebula 

OIL OGLE-II coverage; OVI: OVI i6830 present; P: Other MASH observations; V: FLAMES observations 

References: 

(a) Parker et al. (2006), (b) This work, (c) Jacoby & Van de Steene (2004), (d) Matsunaga, Fukushi & Nakada (2005) 
(e) Zhang & Liu (2003), (f) Gorny (2006), (g) Shaw & Wirth (1985), (h) Sanduleak (1976), (i) Allen (1979) 
(j) Haro (1952), (k) Hajduk et al. (2008), (1) Kohoutek (2002) 

Table 1. PN mimics removed from our PN sample. Columns Rl and R2 are the emission line ratios [OIII] /14363/Hy and [OIII] 
A5001/U/3, respectively (Gutierrez-Moreno, Moreno & Cortes 1995). NIR magnitudes are from 2MASS (Skrutskie et al. 2006). 
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Fig. 6. AAOmega spectrum (left) and OGLE-III time series pho- 
tometry (right) of the newly identified M-dwarf flare star PHR 
1805-2659. 



consistent with intense chromospheric activity. The activity is 
reflected in the OGLE-III photometry as an increasing sequence 
of flares separate to a declining quiescent sequence. 

3.6. Periodic variable determination 

The lomb-scargle technique was used as our main method for 
finding periodic variables (Lomb 1976; Scargle 1982; Press et 
al. 1992). A matlab implementation of the period function (Press 
et al. 1992) was acquired and modified to use pgplot to gener- 
ate periodograms covering a long period window (5-100 d) and 



a short period window (0.1-5 d) under gnu octave0 The peri- 
odograms of all stars within the radius of each PN were visu- 
ally inspected for peaks with a significance level of -1% or less. 
There is a small chance that an unrelated variable star could be 
found within a PN (e.g. Sh 2-71, Frew & Parker 2009), however 
it is relatively uncommon to find more than one variable star per 
PN. We emphasise that if there are any of these objects in our 
search, their assignment as likely or possible CSPN identifica- 
tions means they are treated appropriately in our estimate of the 
binary fraction (see later). We checked the results using both pe- 
riod04 (Lenz & Breger 2004) and less often the IRAF task pdm 
(Stellingwerf 1978) with excellent agreement found amongst all 
three methods. Periods on the order of a hundred days were bet- 
ter searched for visually as discussed in the following section. 

3. 7. Aperiodic variable determination 

We investigated outliers in the magnitude dispersion cr versus 
average magnitude plane for all 'true' CSPN with extracted time 
series photometry. The vast majority of outliers were found to 
be spurious after comparison with surrounding stars. A few ob- 
jects showed real variation but most were symbiotics. Indeed, 
unduly large variations can be used as a good tool in revealing 
symbiotics in conjunction with other data. All photometry was 



http ://w w w. gnu. org/software/octave 
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also inspected visually to identify any unusual variations with 
the time axis 'stretched-out' to better resolve the 'bunched-up' 
Bulge seasons. 

4. Results and Discussion 

4.1. New Binary Central Stars 

Table [2] presents 22 binary CSPN resulting from our searches 
with 19 of these presented here for the first time. We include 
M 3-16, H 2-29 and M 2-19 from Miszalski et al. (2008b). All 
discoveries are from OGLE-III except for JaSt 66 which is an 
OGLE-II discovery. A 'C is placed in the CSPN column of 
Tab. |2] where spectroscopic confirmation of the CSPN has been 
achieved, otherwise the TLP identification is given (see Sec. 
13.21) . We emphasise that the 'T' and 'L' binaries still require 
spectroscopic confirmation. The 'Type' column uses nomencla- 
ture adopted from DM08 to describe the cause of variability as 
either irradiation effects (I), ellipsoidal variation (El) or eclipses 
(Ec). The epoch of primary minimum Eq, the period P and peak- 
peak amplitude A (only for irradiated binaries) were calculated 
using period04 (Lenz & Breger 2004). Periods for the seven 
eclipsing systems were selected to give two minima per cycle as 
usual. When two different minima per cycle were found we put 
this down to ellipsoidal variation. Otherwise we have assumed ir- 
radiation effects are the cause of the variation and select periods 
that show only one minimum per cycle. Future RV observations 
may possibly double the period of a few irradiated variables to 
become ellipsoidal variables, but for now we settle on the peri- 
ods given in Tab. [2 The average OGLE-III /-band magnitudes 
(fhi) are accurate to 0.1-0.2 mag (Udalski et al. 2008). More un- 
certain magnitudes are marked with ':' including K 6-34 where 
we give the average magnitude of the two close stars adding -0.2 
mag of uncertainty (see below) and H 1-33 which is affected by 
moderate nebular contamination. 

Figure [7] depicts images and phased lightcurves for the 21 
periodic binaries from Tab. [2 Ha images measuring 30 x 30 
arcsec 2 are sourced from CCD photometry where available as 
follows. During our ESO visitor mode program 079.D-0764(B) 
we obtained the 30 s Ha exposure of H 1-33 with NTT/EMMI 
on 21 June 2007. Gemini GMOS South (Hook et al. 2004) ac- 
quisition images of Sab 41, M 3-16 and M 2-19 were obtained 
by us using the Ha filter on 24, 25 and 31 July 2008 during 
Program ID GS-2008B-Q-65. Kovacevic & Parker (private com- 
munication) obtained the Ha image of K 6-34 using the Mosaic 
II camera on the CTIO 4-m Blanco telescope on 13 June 2008. 
The K 6-34 Ha image when compared to the OGLE-III V- and 
/-band images shows the bluer southern star of a pair separated 
by 0.9" is the true CSPN. We rule out any physical association 
with the northern star based on radial velocities, but the close 
proximity of these stars means the OGLE-III photometry is in- 
distinguishable between the two (illustrated in Fig. [7] by the red 
lightcurve matched to the fainter star). Ruffle et al. (2004) ac- 
quired Ha images of Bl 3-15 and H 2-29 with NTT/EMMI un- 
der ESO program 67.D-0527(A) that we have smoothed and re- 
produced here. For the remaining 14 PNe we use data from the 
SHS, mostly as quotient images constructed by dividing Ha by 
the contemporaneous broad-band Short-Red exposure (indicated 
by 'Q' in Fig. [7]). For Pe 1-9 we use the Short-Red only because 
the Ha image is saturated. Larger images and a more detailed 
appraisal of nebular morphologies will be presented in another 
paper in this series (Miszalski et al., in prep). 

Figure [8]presents M 1-34 as the only genuine aperiodic vari- 
able found and we consider it a likely binary. Other aperiodic 




Fig. 8. Ha image, OGLE-III /-band image and time series pho- 
tometry of M 1-34. Images are 30 x 30 arcsec 2 with North to top 
and East to left. 



candidates were ruled out as either symbiotic stars (Sec. 13.5.11) 
or were caused by systematic effects (Sec. 13.41) . The Ha im- 
age reproduced here was acquired by Ruffle et al. (2004) with 
NTT/EMMI under ESO program 71.D-0448(A). The star im- 
mediately to the NE of our chosen CSPN is not variable, but we 
cannot rule it out as the CSPN given its slightly more plausible 
central position. Longslit spectroscopy of both stars is required 
to clarify the situation. Further photometric monitoring of this 
object is especially important as the OGLE-III field is no longer 
observed by the survey. 

4.2. Period Distribution and Selection Effects 

Numerous population synthesis calculations have predicted the 
observed period distribution of close binary central stars (de 
Kool 1990, 1992; de Kool & Ritter 1993; Yungelson, Tutukov 
& Livio 1993; Han, Podsiadlowski & Eggleton 1995). Although 
there are many assumptions built into these models, only the ini- 
tial mass ratio distribution (IMRD) and acE markedly influence 
the predicted period distributions. The models therefore offer a 
valuable tool to directly constrain the IMRD and acE when com- 
pared to the observed period distribution of binary CSPN. Until 
now these comparisons have been severely limited by the very 
small sample size and especially the uncertain observational bi- 
ases that constituted the sample (DM08). The well-understood 
observational biases and threefold increase in sample size pro- 
vided by our OGLE-III discoveries addresses both concerns en- 
abling a fresh and more meaningful comparison to be made. 

Figure [9] depicts the period distributions of our new sam- 
ple, the known binaries excluding the uncertain SuWt 2 and Hb 
12 (DM08) and the combined distribution. Their appearance is 
rather sensitive to the bin size for which we have chosen log 
P = 0.1. Table [3] outlines the six different models from de Kool 
(1992) and Han et al. (1995) that have been overlaid on the ob- 
served distribution. Han et al. (1995) models make use of an ad- 
ditional a t h parameter that has the effect of increasing acE when 
a t h = 1.0. Note that the models with lower acE values (dotted 
lines) have steeper slopes past the distribution maximum and that 
the position of the distribution maximum is strongly affected by 
the IMRD, i.e. dK92 M3 and M4 use a random IMRD, while the 
other models use a constant IMRD i.e. dN oc dq (de Kool 1992). 
Other models are available for comparison but these overwhelm- 
ingly utilise a constant IMRD and are fairly well represented by 
some of the models already shown here. 

The models were normalised such that their value equals four 
at log P = -0.5 (similar to DM08) and reduced^ 2 values were 
subsequently calculated to assist comparison with the observed 
distribution (Tab. [3]). An alternative normalisation of the mod- 
els such that the maxima were equal to four was also made, 
but the x 1 results were not sufficiently different to change our 
analysis. With X 2 £ 1, dK92 M3 is the best fit (-80% prob- 
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Fig. 7. Images and phased lightcurves of the 21 confirmed, true and likely periodic binary CSPN. Each row per panel depicts left 
to right: an Ha image of the nebula, an OGLE /-band image, and a phased OGLE /-band lightcurve of the central star (marked on 
the /-band image). Periods and ephemerides are given in Tab. [2] Two indistinguishable lightcurves are given for K 6-34 (see text). 
All images are 30 x 30 arcsec 2 with North to top and East to left. A 'Q' designates an Ha quotient image and 'SR' an Ha off-band 
image (see text). 
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- 


BLG155 


357.6-03.3 


H2-29 




T 


Ec,I? 


3619.52078 


0.244110 


18.1 


- 


BLG101 


000.2-01.9 


M2-19 




C 


El 


4224.78224 


0.670170 


15.9 


- 


BLG172 


358.7-03.0 


K6-34 




C 


El 


3525.59585 


0.393309 


16.4: 


- 


BLG142 


357.0-04.4 


PHR 1756- 


-3342 


C 


I 


4331.41964 


0.265733 


18.0 


551 


BLG214 


001.8-02.0 


PHR 1757- 


-2824 


T 


Ec,I 


4286.62547 


0.799209 


18.1 


453 


BLG104 


001.2-02.6 


PHR 1759- 


-2915 


L 


Ec 


4524.84332 


1.103664 


17.5 


- 


BLG188 


000.5-03. la 


MPA 1759 


-3007 


L 


El 


3930.67937 


0.503604 


18.2 


- 


BLG215 


001.9-02.5 


PPA 1759- 


2834 


L 


I 


3099.73716 


0.305848 


18.2 


501 


BLG135 


357.1-05.3 


BMP 1800 


-3408 


T 


EcJ 


2934.58226 


0.144777 


18.1 


80 


BLG189 


000.9-03.3 


PHR 1801- 


-2947 


L 


I 


2132.58866 


0.315998 


17.5 


110 


BLG151 


357.9-05.1 


M 1-34 




L 


Ec 






16.8 




BLG233 


003.1-02.1 


PHR 1801- 


-2718 


L 


I 


4339.47776 


0.321972 


17.3 


283 


BLG196 


001.8-03.7 


PHR 1804- 


-2913 


C 


I 


3578.83476 


6.659968 


15.0 


13 


BLG241 


004.0-02.6 


PHR 1804- 


-2645 


T 


Ec 


3419.84190 


0.624505 


18.5 





Table 2. The 22 confirmed, true and likely binary CSPN discovered in this work and Miszalski et al. (2008b). See text for column 
descriptions. 




log P 



Fig. 9. Orbital period distribution of new binaries (shaded, this 
work), known binaries (hatched, DM08) and the total sample 
(thick line). Details of the overlaid models and their normalisa- 
tion may be found in Tab. [3] and in the text. 



ability), whereas dK92 M2 (-45%) and dK92 M4 (-22%) are 
also acceptable despite the nonaligned maxima between the ob- 
served and model distributions. The disagreement is significant 
for dK92 Ml and H95 S4 (-2%) and highly significant for H95 
S6 (<0.05%). 

Realising dK92 M3 as the best fit model strengthens con- 
siderably earlier suspicions that prefer a random IMRD based 
on a much smaller sample than that used here (de Kool 1990; 
de Kool & Ritter 2003). A random IMRD is supported by the 
observed predominance of late-type secondaries (DM08) that is 
predicted by equivalent models (see Fig. 10c, de Kool & Ritter 
1993; Politano & Weiler 2006) and the secondary mass distribu- 
tion of spectroscopic binaries studied by Goldberg et al. (2003). 



Model 


IMRD 


cuce 




x 1 


dK92 Ml 


dN oc dq 


1.0 




1.65 


dK92 M2 


dN oc dq 


0.3 




1.02 


dK92 M3 


M 2 random from IMF 


1.0 




0.72 


dK92 M4 


M 2 random from IMF 


0.3 




1.17 


H95 S4 


dN oc dq 


1.0 


0.0 


1.66 


H95 S6 


dN oc dq 


1.0 


1.0 


2.54 



Table 3. List of models used in comparison to the observed 
distribution with their adopted initial mass ratio distributions 
(IMRD), acE and a t u values. The result of a reduced^ 2 fit to the 
observed period distribution is also given. Models prefixed with 
dK92 and H95 are from de Kool (1992) and Han et al. (1995), 
respectively. 



We refrain from drawing any strong conclusions based on our 
comparison, but nevertheless it is concerning that the majority 
of literature models assume a constant IMRD that predicts sub- 
stantially more binaries at longer periods. A reliance on these 
models has created uncertainty regarding whether the large dis- 
crepancy for log P ^ is intrinsically real or caused by selection 
effects (DM08). DM08 compared the then known period distri- 
bution with the H95 S4 and S6 models to suggest a severe inade- 
quacy in all CE population synthesis models. Our larger sample 
and usage of the best fit dK92 M3 reduces the severe discrep- 
ancy to only a modest one. The more modest discrepancy is now 
more easily explained by selection effects rather than a severe 
flaw in our understanding of CE evolution suggested by DM08. 

At very short periods (log P ^ -1) the models are in good 
agreement with observations and any deficiency could be at- 
tributed to the sampling frequency of each field (see Sec. 14.3b . 
In assessing the contribution of selection effects towards the 
long period discrepancy (log P ^ 0), a very useful compari- 
son can be made with the observed period distribution of post- 
CE white dwarf/main sequence (WDMS) binaries (Rebassa- 
Mansergas et al. 2008, hereafter RM08). In spite of the small 
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samples and differences in selection effects associated with the 
discovery method of post-CE WDMS (spectroscopic and photo- 
metric, Schreiber & Gansicke 2003) and CSPN (photometric) bi- 
naries, the distributions bear a very strong resemblance to each 
other which strongly suggests they are an acceptable depiction 
of the true post-CE binary population. RM08 reach the same 
conclusion based on Monte Carlo simulations that showed they 
should have been sensitive to detecting longer period binaries 
in their search. Similarly, DM08 claimed the then known CSPN 
binaries were also representative of the true population, despite 
serious misgivings about the contribution of selection effects to 
the sample discovery. DM08 supported their claim by conduct- 
ing a theoretical investigation into the expected amplitudes of 
irradiated binaries for periods P ^ 3 days and concluded many 
systems should have been detected by BOO down to amplitudes 
of -0.1 mag. 

If the equivalence of the CSPN and post-CE WDMS period 
distributions is indeed indicative of the true post-CE popula- 
tion, then selection effects alone may not account for the mod- 
est discrepancy at longer periods and a revision of the popula- 
tion synthesis models may be required. Alternatively, the most 
likely cause of the discrepancy is the / ~ 20 magnitude limit 
of OGLE-III. At 8 kpc we are likely insensitive to the detec- 
tion of M-dwarf secondaries at / ^ 21 (Kirkpatrick & McCarthy 
1994), but irradiation effects may permit detection of some sys- 
tems. Some secondaries are expected to have spectral types later 
than M (de Kool & Ritter 1993; Politano & Weiler 2006), but 
these are likely to be rare (Farihi, Becklin & Zuckerman 2005). 
A dominant bright primary may also preclude faint secondaries 
from being detected in some cases. Some more non-detections 
could be non-irradiated eclipsing binaries, i.e. if M 3-16 were 
observed at low-moderate inclinations, however these binaries 
would not necessarily be restricted to log P ^ 0. A stronger test 
of the discrepancy could come from RV monitoring surveys that 
target the longer period regime up to tens of days (Mendez 1989; 
De Marco et al. 2004; Sorenson & Pollacco 2004; Afsar & Bond 
2005), but the results from these surveys are largely inconclusive 
with few, if any, firm periods found as the target sampling seems 
far from ideal and the typical number of 10-20 observations per 
object is still quite small. Indeed, the high levels of variability 
seen (excluding Mendez 1989) could be attributed to winds and 
pulsations from the central star alone, rather than being caused 
by binary companion. Surveys with more frequent sampling may 
be able to fill the modest discrepancy at longer periods, but we 
suspect the number of binaries with periods of a few days would 
remain small based on the post-CE WDMS distribution which 
includes binaries found using this method. 

The inital investigation into expected variability amplitudes 
by DM08 raises perhaps the most important selection effect of 
sensitivity. Around 30% of OGLE-III binaries exhibit ampli- 
tudes ^0.1 mag which is markedly different from the known 
binaries where only Sp 1 has such an amplitude. Although we 
have considerably improved sensitivity compared to BOO in the 
low amplitude regime, the fact that we still find a large absence 
of binaries at longer periods reinforces the representative nature 
of the current period distribution. Alternatively, some low am- 
plitude binaries may be lost in the larger intrinsic scatter in some 
fields (Sec. 13.4b . Further investigations following DM08 into the 
expected low-amplitude population for a variety of irradiated bi- 
nary configurations and inclinations would be beneficial. 

Constraining acE is more problematic than constraining the 
IMRD. Even though acE = 1.0 was used in calculating the 
best fitting dK92 M3, we caution against adopting acE = 1.0 
from this work, mainly because different definitions of acE in 



the literature complicate comparisons between different models 
(Livio 1996) and suspicions that acE may not be constant (Livio 
1996; Politano & Weiler 2007). Furthermore, measuring a post- 
maximum gradient from the observed period distribution is still 
hindered by the small sample size even with the new OGLE dis- 
coveries. 

4.3. Binary Fraction 

We are now able to derive the first independent estimate of the 
short period binary CSPN fraction since BOO. However, not all 
OGLE-III fields are sampled uniformly and this can unduly bias 
any estimate of the binary fraction. Understanding the capacity 
of each field to detect short binary CSPN periods (P ^ 1 d) is 
required before any less sensitive fields can be identified and re- 
moved. Usually the Nyquist frequency would be used but this 
is not necessarily defined for each field that has been sampled 
intermittently over the course of each Bulge season and many 
years. We therefore calculated the number of consecutive obser- 
vations per field, N(At\ ) and Af(Af me d), for separation intervals 
At\ < 0.5 d and 0.5 d < A^ me( j < 1.5 d, respectively. These inter- 
vals were chosen because their Nyquist frequencies lied within 
our desired detection range. Sensitivity to longer periods up to 
hundreds of days is provided by the long time baseline of the 
data and because of limited sensitivity to P ^ 0.1 d our search 
was only conducted for longer periods. 

Figure[lOldisplays the relation between N(At\ ) and N(At mQ ^) 
for all OGLE-III fields. Note the excellent sampling overall and 
that the new binaries lie only in well- sampled fields along the 
main locus. The data suggest fields located within the rather con- 
servatively chosen shaded region be excluded from binary frac- 
tion calculations (Fig. [TO]). This region contains the field that ob- 
served BMP 1800-3408 which we consider to be an anomalous 
detection. Figure [TT] shows the spatial distribution of N(At me ^) 
amongst OGLE-III fields with new binaries overlaid. Note the 
strong concentration towards better sampled southern latitude 
fields with Pe 1-9 being the only binary detected at positive lati- 
tudes. 

To estimate the binary fraction we make cuts to Tab. IA.1I 
by first removing the PNe belonging to fields largely insensitive 
to binary CSPN detection. We further remove PNe with remarks 
'NULL', 'NEB' or 'P' as we are not able to examine the variabil- 
ity or otherwise of these objects (Sec. 13.21) . These removed ob- 
jects contribute towards a 'rejected' sample that may include ob- 
jects we suspect of being affected by systematic effects labelled 
'S(?)\ This leaves two samples, namely the smaller sample con- 
taining only 'true' CSPN identifications (T) and a larger sample 
containing 'true' and 'likely' CSPN identifications (T+L). Some 
binaries located in fainter PNe fall into the 'likely' category in 
line with other faint PNe and are not included in the T sample. 
The T+L sample constitutes a lower limit to the binary fraction. 
Note also that we exclude M 1-34 given its quite different nature 
to the short period binaries. The binary fraction is calculated as 
the number of binaries (BIN) present in each sample with respect 
to the total number of objects that were not rejected. 

Table @] contains the calculated binary fractions for both T 
and T+L samples with and without the 10 PNe affected by sys- 
tematic effects. Our close binary fraction of 12-21% is in good 
agreement with 10-15% found by BOO and 12-33% determined 
by Frew & Parker (2009). However, our sample sizes of ~65 PNe 
(T) and -145 (T+L) PNe only allow for a meaningful compar- 
ison with the BOO sample of -100 PNe. The early results from 
Frew & Parker (2009) are based on a sample of only 33 PNe se- 
lected from a larger catalogue of nearby PNe (Frew 2008) which 
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Fig. 10. The relation between the number of consecutive obser- 
vations N(At\ ) and N(At mQd ) of OGLE-III fields, where At\ < 
0.5 d and 0.5 d < At mQ a < 1.5 d, respectively. Fields that contain 
the newly discovered OGLE-III binaries from Fig. [7] (open cir- 
cles) are located within a locus of well- sampled fields (top). The 
shaded region at lower left contains fields excluded from binary 
fraction calculations. 
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Galactic Longitude 

Fig. 11. Same as Fig. Q] but with PNe replaced by binaries from 
Tab. [2] (circled points) and colours indicative of the number of 
consecutive observations N(At mQ d) separated by the interval 0.5 
d < A^ m ed < 1.5 d. Less well-sampled fields coloured grey are 
excluded from binary fraction calculations. 



has considerable overlap with the BOO sample. Our estimate also 
agrees well with model predictions (e.g. -17%, Han et al. 1995) 
and with the fraction of post-CE WDMS binaries (e.g. -15%, 
RM08). Accounting for our insensitivity to M-type secondaries 
(Sec. 14.21) could lead to an upward revision of our estimate to 
-35%, but this is left for more sensitive future surveys to con- 
firm. If we include fields removed on the basis of our sampling 
cut, then this slightly reduces our derived binary fraction to - 10- 
18%. The large range in our estimate is caused by the lack of 
OGLE-III V-band data that would have otherwise enabled more 
'likely' identifications to become 'true' via appropriate colour 
selections. However, we note the lower limit is considered quite 
robust given that all CSPN candidates for the 'likely' sample 
were examined for variability. 

5. Conclusion 

A large catalogue of -300 Galactic Bulge PNe was constructed 
after removing 24 new symbiotic stars with the assistance of 



Sample 


S 


Passed 


Rejected 


Binaries 


Fraction (%) 


T 


N 


61 


149 


13 


21.3 


T+L 


N 


139 


149 


18 


12.9 


T 


Y 


71 


139 


13 


18.3 


T+L 


Y 


149 


139 


18 


12.1 



Table 4. Binary fraction calculations for the two samples. The 
'S' column indicates whether PNe affected by systematic effects 
are included. See text for further details. 



deep optical spectroscopy. Time series photometry of identified 
central stars were extracted from the OGLE-III survey and anal- 
ysed for periodic and aperiodic variability. Nebular contamina- 
tion and other systematic effects were discussed in the context of 
often non-stellar PNe in a stellar-optimised survey. Images and 
lightcurves of 21 new periodic close binaries and 1 new likely 
aperiodic binary were presented, more than doubling the num- 
ber of known close binary CSPN. 

The observed period distribution was compared with pre- 
dictions made by CE population synthesis models. A best fit 
was achieved for model 3 from de Kool (1992), which adopted 
acE = 1.0 and a random IMRD. Only a modest discrepancy be- 
tween the best fit model and the observed distribution is found 
at longer periods, contrary to most other models that produce a 
more severe discrepancy. Further comparison with the post-CE 
WDMS binary period distribution suggests the CSPN and post- 
CE WDMS binary period distributions are an adequate depiction 
of the real post-CE population (i.e. binaries with periods longer 
than about one day are rare). If the discrepancy cannot be re- 
solved by photometric surveys more sensitive to M- and L-type 
secondaries or RV surveys of large numbers of CSPN conducted 
with more frequent sampling, then a revision of the CE popula- 
tion synthesis models may be required. 

We have firmly established the close binary fraction to be 
10-20% from our larger and entirely independent sample of 
binary PNe. This agrees well with previous photometric (10- 
15%, BOO) and CE population synthesis models (17%, Han et 
al. 1995). RV monitoring programs have proved rather unreli- 
able and inconclusive, however our results agree well with an 
estimate of 15% by Mendez (1989). Our estimate does not sup- 
port the hypothesis that most PNe derive from binaries (Moe 
& De Marco 2006), nor the high levels of variability mistaken 
for binarity in RV monitoring surveys (De Marco et al. 2004; 
Sorenson & Pollacco 2004; Afsar & Bond 2005). With the close 
binary fraction now well determined, more effort is required at 
the longer period end through detailed study of CSPN showing 
composite spectra (e.g. Sorenson & Pollacco 2004), resolved bi- 
naries (Ciardullo et al. 1999) and infra-red excesses (e.g. Frew 
& Parker 2009). As a shaping mechanism of nebular morpholo- 
gies close binaries seem to play an important but not major role. 
Indeed, it has been suggested that the largest role may be ful- 
filled by low-mass or sub-stellar secondaries (Soker 1997) and 
the aforementioned studies on longer period binaries will help 
in this area. This important topic will be further discussed in a 
following paper with an examination of nebular morphologies of 
our new sample. 

Emphasis should now be placed more on radial velocity 
monitoring programs and an AAOmega survey of brighter CSPN 
covered by OGLE-III would prove particularly useful in better 
understanding the selection effects of our survey. However, there 
are still advancements to be made via photometric campaigns. 
A bias against the brightest PNe in OGLE-III requires dedic- 
tated independent followup perhaps using narrow-band filters to 
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remove nebular contamination. Planned surveys such as Vista 
Variables in the Via Lactea (VVV, see Arnaboldi et al. 2007) are 
likely to find more binaries outside OGLE-III fields provided 
field sampling is not too heavily biased towards higher stellar 
density fields already covered by OGLE-III. Since VVV is con- 
ducted in the A^-band the survey may be more sensitive to bina- 
ries with later-type secondaries than OGLE-III and would there- 
fore provide greater understanding of the selection effects of our 
survey. Dedicated photometric surveys with telescopes such as 
Sky Mapper (Keller et al. 2007) are required to survey regions 
outside microlensing surveys or where the sampling is not al- 
ways ideal and such surveys would prove very fruitful indeed in 
monitoring the now very large population of Galactic Bulge and 
Southern Galactic Plane PNe (Parker et al. 2006; Miszalski et al. 
2008a). 
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Table A.l. Planetary nebulae with OGLE-III coverage. Remarks are explained in the text. 
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T 
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4 
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17 


51 
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39 


13 
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17 


51 


23, 


,0 


-35 


31 
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1 
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17 


51 
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,9 
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59 


54 


L 
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.3+02, 


.2 
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17 


51 
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o 
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23 


18 


L 
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.8 
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17 


51 
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.2 
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12 


47 


p 
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.5 
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17 


51 
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.6 
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36 


01 


T 
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,7 
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17 


51 
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,7 
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02 


34 


T 
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.6 
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17 


51 
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,6 
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47 


36 


L 


BLG171 
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4 
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17 


52 
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.2 
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17 


50 


T 
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.1 
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17 


52 
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,7 
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44 


08 


L 
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.3 
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17 


52 
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,8 
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16 


42 
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.8 
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17 


52 
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.0 
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05 


14 


T 
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.0 


PHR 1752- 
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17 


52 
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.8 
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33 


08 


T 
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.6- 
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.Oa 
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17 


52 
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,1 
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30 


07 


TB 
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355, 
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.3 
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17 


52 
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.0 


-35 


04 


18 


L 
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.6 
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17 


52 
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.2 
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30 
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T 
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.4- 
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.2 


M 2-16 




17 


52 
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.3 
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45 


51 
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BLG195 
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.6- 
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.3 
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17 


52 
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.0 
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06 


39 
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358, 
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.5 
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17 


52 


36, 


.5 
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16 


27 
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006 


.4+02, 


.0 
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17 


52 
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,5 


-22 


21 


57 


NULL 


BLG101 


000, 


.0- 
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.8 


JaSt 83 




17 


52 
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1 
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50 


59 


L 


BLG180 
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.3 
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17 


52 
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49 


35 
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52 
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Field 


PNG 




Name 




RA (J2000) 


DEC (J2000) 


Remarks 


BLG133 


356.5- 


-04.1 


PPA 1754- 


■3358 


17 


54 


03.1 


-33 


58 


51 


NULL 


BLG101 


OOO.O- 


-02.1 


MPA 1754 


-2957 


17 


54 


04.3 


-29 


57 


27 


L 


BLG119 


355.5- 


-04.8 


PHR 1754- 


-3515 


17 


54 


17.5 


-35 


15 


39 


L 


BLG101 


000.4- 


-01.9 


M2-20 




17 


54 


25.4 


-29 


36 


09 


T,S? 


BLG119 


355.2- 


-05.0 


PHR 1754- 


-3533 


17 


54 


30.7 


-35 


33 


08 


L 


BLG133 


356.2- 


-04.4 


Cn 2-1 




17 


54 


32.9 


-34 


22 


22 


NULL 


BLG172 


358.7- 


-03.0 


K6-34 




17 


54 


41.3 


-31 


31 


43 


T,B,BIN 


BLG172 


359.1- 


-02.9 


M3-46 




17 


55 


05.5 


-31 


12 


17 


NEB 


BLG172 


359.3- 


-02.8 


MPA 1755 


-3058 


17 


55 


13.1 


-30 


58 


13 


L 


BLG103 


000.2- 


-02.3 


Bl 3-10 




17 


55 


20.6 


-29 


57 


36 


NEB 


BLG181 


359.7- 


-02.6 


H 1-40 




17 


55 


36.0 


-30 


33 


33 


T 


BLG103 


000.0- 


-02.5 


K6-36 




17 


55 


52.8 


-30 


15 


41 


T 


BLG164 


358.8- 


-03.3 


PK 358-033 


17 


56 


02.4 


-31 


38 


23 


T 


BLG205 


000.9- 


-02.0 


Bl 3-13 




17 


56 


02.6 


-29 


11 


16 


NULL 


BLG103 


359.9- 


-02.6 


PHR 1756- 


-3019 


17 


56 


06.5 


-30 


19 


36 


P 


BLG156 


357.9- 


-03.8 


H2-30 




17 


56 


13.9 


-32 


37 


22 


T 


BLG172 


359.3- 


-03.1 


M 3-17 




17 


56 


25.7 


-31 


04 


17 


T,S? 


BLG172 


359.2- 


-03.1 


PHR 1756- 


-3112 


17 


56 


28.8 


-31 


12 


35 


L 


BLG214 


001.5- 


-01.8 


JaSt 2-19 




17 


56 


33.8 


-28 


30 


30 


L 


BLG205 


001.2- 


-02.0 


PHR 1756- 


-2857 


17 


56 


36.5 


-28 


57 


18 


P 


BLG142 


357.0- 


-04.4 


PHR 1756- 


-3342 


17 


56 


39.6 


-33 


42 


31 


T,B,BIN 


BLG164 


358.6- 


-03.6 


PHR 1756- 


-3157 


17 


56 


43.7 


-31 


57 


33 


P 


BLG134 


356.6- 


-04.7 


PHR 1756- 


-3414 


17 


56 


48.2 


-34 


14 


32 


L 


BLG134 


356.7- 


-04.7 


MPA 1757- 


-3410 


17 


57 


00.7 


-34 


10 


41 


T,B 


BLG103,182 


000.0- 


-02.9 


MPA 1757- 


-3021 


17 


57 


14.2 


-30 


21 


53 


L 


BLG134 


356.7- 


-04.8 


H 1-41 




17 


57 


19.0 


-34 


09 


49 


T 


BLG142 


357.2- 


-04.5 


H 1-42 




17 


57 


25.2 


-33 


35 


43 


T 


BLG214 


001.8- 


-02.0 


PHR 1757- 


-2824 


17 


57 


42.2 


-28 


24 


07 


T,B,BIN 


BLG188 


000.3- 


-02.8 


M 3-47 




17 


57 


43.4 


-30 


02 


30 


T 


BLG165 


358.9- 


-03.6 


PPA 1757- 


3144 


17 


57 


48.7 


-31 


44 


05 


P 


BLG173 


359.4- 


-03.3 


PHR 1757- 


-3106 


17 


57 


51.6 


-31 


06 


11 


L 


BLG104 


000.7- 


-02.7 


M2-21 




17 


58 


09.6 


-29 


44 


20 


NULL 


BLG165 


358.9- 


-03.7 


H 1-44 




17 


58 


10.6 


-31 


42 


56 


T 


BLG173 


359.4- 


-03.4 


H2-33 




17 


58 


12.5 


-31 


08 


06 


T 


BLG143 


357.1- 


-04.7 


H 1-43 




17 


58 


14.4 


-33 


47 


37 


T 


BLG205 


001.6- 


-02.2 


MPA 1758 


-2841 


17 


58 


14.8 


-28 


41 


30 


P 


BLG188 


000.4- 


-02.9 


M3-19 




17 


58 


19.4 


-30 


00 


40 


T 


BLG150 


357.8- 


-04.4 


PHR 1758- 


-3304 


17 


58 


25.9 


-33 


04 


59 


L 


BLG223 


002.6- 


-01.7 


PHR 1758- 


-2729 


17 


58 


28.9 


-27 


29 


40 


P 


BLG223 


002.5- 


-01.7 


Pe 2-11 




17 


58 


31.2 


-27 


37 


06 


T 


BLG143 


357.4- 


-04.6 


M2-22 




17 


58 


32.6 


-33 


28 


37 


T 


BLG223 


002.3- 


-01.9 


PHR 1758- 


-2756 


17 


58 


35.0 


-27 


56 


56 


P 


BLG157,158 


358.5- 


-04.2 


H 1-46 




17 


59 


02.4 


-32 


21 


44 


NULL 


BLG104 


001.2- 


-02.6 


PHR 1759- 


-2915 


17 


59 


02.9 


-29 


15 


01 


L,BIN 


BLG188 


000.5- 


-03.1 


KFL 1 




17 


59 


15.6 


-30 


02 


48 


L 


BLG215 


002.1- 


-02.2 


M3-20 




17 


59 


19.4 


-28 


13 


48 


T 


BLG173,174 


359.5- 


-03.7 


MPA 1759 


-3116 


17 


59 


25.5 


-31 


16 


54 


P 


BLG188 


000.5- 


-03.1a 


MPA 1759 


-3007 


17 


59 


25.6 


-30 


07 


15 


L,BIN 


BLG232 


003.1- 


-01.6 


PHR 1759- 


-2706 


17 


59 


26.2 


-27 


06 


34 


NULL 


BLG232 


003.0- 


-01.7 


PHR 1759- 


-2712 


17 


59 


33.1 


-27 


12 


50 


T,B 


BLG143 


357.7- 


-04.8 


BMP 1759 


-3321 


17 


59 


45.2 


-33 


21 


13 


T,B 


BLG183 


000.2- 


-03.4 


PHR 1759- 


-3030 


17 


59 


47.8 


-30 


30 


35 


P 


BLG215 


001.9- 


-02.5 


PPA 1759- 


-2834 


17 


59 


52.6 


-28 


34 


47 


L,BIN 


BLG151 


358.0- 


-04.6 


Sa 3-107 




17 


59 


55.0 


-32 


59 


12 


T 


BLG232 


003.3- 


-01.6 


PHR 1759- 


-2651 


17 


59 


55.2 


-26 


51 


49 


L 


BLG232 


002.9- 


-01.8 


MPA 1759 


-2719 


17 


59 


55.5 


-27 


19 


17 


P 


BLG166 


359.0- 


-04.1 


M3-48 




17 


59 


56.6 


-31 


54 


27 


NEB 



B. Miszalski et al.: Binary Planetary Nebulae Nuclei towards the Galactic Bulge 



Field 


PNG 




Name 




RA (J2000) 


DEC (J2000) 


Remarks 


BLG188 


000.6- 


-03.2 


MPA 1759 


-3004 


17 


59 


56.7 


-30 


04 


28 


L 


BLG206 


001.6- 


-02.6 


PHR 1759- 


-2853 


17 


59 


58.8 


-28 


53 


57 


P 


BLG158 


358.4- 


-04.5 


K6-38 




17 


59 


59.0 


-32 


35 


58 


T 


BLG183 


000.2- 


-03.4a 


MPA 1800 


-3026 


18 


00 


00.0 


-30 


26 


08 


P 


BLG206 


001.7- 


-02.6 


PPA 1800- 


-2846 


18 


00 


00.7 


-28 


46 


27 


NULL 


BLG183 


000.3- 


-03.4 


MPA 1800 


-3023 


18 


00 


11.1 


-30 


23 


49 


T,B 


BLG135 


356.8- 


-05.4 


H2-35 




18 


00 


18.2 


-34 


27 


40 


L 


BLG215 


002.1- 


-02.4 


PPA 1800- 


■2818 


18 


00 


18.7 


-28 


18 


35 


L 


BLG215 


002.0- 


-02.5 


PPA 1800- 


■2826 


18 


00 


18.7 


-28 


26 


08 


L 


BLG206 


001.5- 


-02.8 


PPA 1800- 


■2904 


18 


00 


22.3 


-29 


04 


39 


T 


BLG135 


357.1- 


-05.3 


BMP 1800 


-3408 


18 


00 


26.1 


-34 


08 


03 


T,B,BIN 


BLG206 


001.8- 


-02.7 


PHR 1800- 


-2842 


18 


00 


36.0 


-28 


42 


51 


P 


BLG105 


001.2- 


-03.0 


H 1-47 




18 


00 


37.7 


-29 


21 


51 


NULL 


BLG232 


003.4- 


-01.8 


PHR 1800- 


-2653 


18 


00 


42.2 


-26 


53 


37 


L 


BLG144 


357.5- 


-05.1 


PPA 1800- 


■3341 


18 


00 


43.2 


-33 


41 


59 


T 


BLG215 


002.2- 


-02.5 


KFL2 




18 


00 


59.8 


-28 


16 


19 


L 


BLG224 


002.8- 


-02.2 


Pe 2-12 




18 


01 


10.1 


-27 


38 


26 


NEB 


BLG189 


000.9- 


-03.3 


PHR 1801- 


-2947 


18 


01 


13.4 


-29 


47 


00 


L,BIN 


BLG151 


357.9- 


-05.1 


M 1-34 




18 


01 


22.1 


-33 


17 


43 


L,BIN? 


BLG233 


003.1- 


-02.1 


PHR 1801- 


-2718 


18 


01 


24.5 


-27 


18 


09 


L,BIN 


BLG224 


002.7- 


-02.4 


PPA 1801- 


■2746 


18 


01 


32.4 


-27 


46 


07 


L 


BLG216 


002.1- 


-02.8 


PHR 1801- 


-2831 


18 


01 


32.4 


-28 


31 


45 


L 


BLG216 


002.2- 


-02.7 


M2-23 




18 


01 


42.7 


-28 


25 


44 


T,S 


BLG151 


358.0- 


-05.1 


Pe 1-11 




18 


01 


42.7 


-33 


15 


25 


NULL 


BLG216 


002.4- 


-02.6 


PHR 1801- 


-2809 


18 


01 


45.4 


-28 


09 


37 


P 


BLG224 


002.6- 


-02.5 


MPA 1801 


-2755 


18 


01 


50.1 


-27 


55 


26 


P 


BLG216 


002.5- 


-02.6 


MPA 1802 


-2803 


18 


02 


04.5 


-28 


03 


39 


L 


BLG189 


000.7- 


-03.7 


M3-22 




18 


02 


19.2 


-30 


14 


25 


T 


BLG144 


357.5- 


-05.5 


PPA 1802- 


-3350 


18 


02 


23.8 


-33 


50 


48 


L 


BLG105 


001.4- 


-03.4 


ShWi 1 




18 


02 


26.4 


-29 


25 


05 


L 


BLG175 


359.7- 


-04.4a 


PPA 1802- 


■3124 


18 


02 


30.3 


-31 


24 


17 


P 


BLG159 


359.0- 


-04.8 


M2-25 




18 


02 


46.6 


-32 


09 


29 


P 


BLG207 


002.0- 


-03.1 


PHR 1802- 


-2847 


18 


02 


48.5 


-28 


47 


41 


L 


BLG207 


001.8- 


-03.2 


MPA 1802 


-2900 


18 


02 


50.2 


-29 


00 


20 


L 


BLG224 


003.0- 


-02.6 


KFL4 




18 


02 


51.4 


-27 


41 


01 


T 


BLG216 


002.6- 


-02.8a 


MPA 1802 


-2807 


18 


02 


53.6 


-28 


07 


56 


L 


BLG216 


002.6- 


-02.8 


PHR 1803- 


-2804 


18 


03 


08.4 


-28 


04 


15 


P 


BLG233 


003.5- 


-02.3 


PHR 1803- 


-2702 


18 


03 


09.8 


-27 


02 


46 


L 


BLG207 


002.0- 


-03.2 


PHR 1803- 


-2848 


18 


03 


09.8 


-28 


48 


51 


L 


BLG233 


003.6- 


-02.3 


M2-26 




18 


03 


12.0 


-26 


58 


30 


L 


BLG159 


359.0- 


-04.9 


PHR 1803- 


-3218 


18 


03 


16.8 


-32 


18 


26 


L 


BLG184 


000.3- 


-04.2 


MPA 1803 


-3043 


18 


03 


21.7 


-30 


43 


36 


P 


BLG216 


002.4- 


-03.1 


PPA 1803- 


■2826 


18 


03 


24.7 


-28 


26 


25 


NULL 


BLG184 


000.3- 


-04.2a 


MPA 1803 


-3046 


18 


03 


25.8 


-30 


46 


25 


P 


BLG159 


358.7- 


-05.1 


SB 53 




18 


03 


28.6 


-32 


37 


25 


L 


BLG190 


001.2- 


-03.8 


PHR 1803- 


-2947 


18 


03 


31.2 


-29 


47 


11 


L 


BLG207 


002.0- 


-03.4 


PPA 1803- 


■2855 


18 


03 


37.4 


-28 


55 


41 


L 


BLG175 


359.9- 


-04.5 


M2-27 




18 


03 


52.6 


-31 


17 


47 


T,S 


BLG159 


358.7- 


-05.2 


H 1-50 




18 


03 


53.5 


-32 


41 


42 


T,S? 


BLG234 


003.8- 


-02.4 


PHR 1804- 


-2653 


18 


04 


02.6 


-26 


53 


27 


L 


BLG136 


357.1- 


-06.1 


M3-50 




18 


04 


05.3 


-34 


28 


38 


NEB 


BLG216,217 


002.4- 


-03.2 


Sa 3-115 




18 


04 


05.5 


-28 


27 


51 


L 


BLG167 


359.6- 


-04.8 


H2-36 




18 


04 


07.7 


-31 


39 


10 


L 


BLG196 


001.6- 


-03.7 


MPA 1804 


-2926 


18 


04 


10.9 


-29 


26 


33 


P 


BLG217 


002.6- 


-03.1 


PHR 1804- 


-2816 


18 


04 


13.2 


-28 


16 


12 


P 


BLG152 


358.4- 


-05.5 


PHR 1804- 


-3306 


18 


04 


14.9 


-33 


06 


30 


T,B 


BLG225 


002.9- 


-03.0 


PHR 1804- 


-2757 


18 


04 


17.3 


-27 


57 


12 


L 
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Field PNG Name RA (J2000) DEC (J2000) Remarks 



BLG196 

x-/ i j\ — j i y \j 


001. 


.8- 


-03. 


j 


PHR 1 804- 

x xxxv i u vy r 


-2913 


18 


04 


28 


.5 


-29 


13 


57 


T,BIN 


BLG208 


002. 


.3- 


-03, 


A 


H 2-37 




18 


04 


28 


.8 


-28 


37 


38 


T 


BLG127 


356. 


,7- 


-06. 


,4 


H 1-51 




18 


04 


29, 


.3 


-34 


58 


00 


L 


BLG175 


000. 


.2- 


-04. 


.6 


Sa 3-117 




18 


04 


44 


.2 


-31 


02 


49 


p 


BLG190 


000. 


.9- 


-04. 


.2 


PHR 1804- 


-3016 


18 


04 


48, 


.0 


-30 


16 


49 


L 


BLG217 

X-J 1 J > — J ^ 1 / 


002. 


.4- 


-03. 


,4 


PHR 1 804- 

x xxxv i yj \J r 


-2833 


18 


04 


55. 


.4 


-28 


33 


08 


PB 


BLG153 


358. 


.6- 


-05. 


.5 


M 3-51 




18 


04 


56, 


.2 


-32 


54 


01 


L 


BLG241 


004. 


0- 


-02. 


.5 


PHR 1 804- 

X XXXV 1 KJ u~ 


-2642 


18 


04 


58. 


.8 


-26 


42 


35 


p 


BLG241 

X-J 1 > — J ^ T 1 


004. 


.0- 


-02. 


.6 


PHR 1 804- 

X XXXV 1 UUT 


-2645 


18 


04 


59, 


.5 


-26 


45 


17 


T,BIN 


BLG176 


000. 


.3- 


-04. 


.6 


M 2-28 




18 


05 


02, 


.6 


-30 


58 


17 


NEB 


BLG234 


003. 


.6- 


-02. 


.8 


MPA 1805 

_LVXX xX X kjKJ^J 


-2712 


18 


05 


03, 


.2 


-27 


12 


33 


p 


BLG185 

X-J 1 > — J X 


000. 


.6- 


-04. 


.5 


PM 1 -206 

x ivi x ^ \j \j 




18 


05 


12 


.5 


-30 


42 


13 


TB 


BLG234 


003. 


.5- 


-02. 


.9 


MPA 1805 

1V11 1 X X UV7 J 


-2721 


18 


05 


13 


1 


-27 


21 


08 


NULL 


BLG241 


004. 


.2- 


-02. 


.5 


PHR 1805- 


-2631 


18 


05 


20, 


.2 


-26 


31 


45 


L 


BLG217 


002. 


.6- 


-03. 


.4 


M 1-37 




18 


05 


25, 


,7 


-28 


22 


04 


T 


BLG234 


004. 


.0- 


-02, 


.7 


PPA 1805- 

X X xX j^yjyj^J 


■2649 


18 


05 


26 


.4 


-26 


49 


03 


p 


BLG168 


359. 


,7- 


-05.0 


PHR 1805- 


3140 


18 


05 


30, 


.2 


-31 


40 


16 


L 


BLG217 


002. 


.9- 


-03. 


.3 


PHR 1805- 

x xxxv i yj \J 


-2804 


18 


05 


44, 


.4 


-28 


04 


46 


L 


BLG208 


002. 


1- 


-03. 


.8 


MPA 1805 


-2902 


18 


05 


55, 


,4 


-29 


02 


46 


p 


BLG241 


004. 


.3- 


-02. 


.6 


H 1-53 




18 


05 


57, 


.4 


-26 


29 


42 


NULL 


BLG234 


003. 


.6- 


-03. 


.0 


PHR 1805- 

X XXXV 1UUJ 


-2723 


18 


05 


57, 


.8 


-27 


23 


09 


L 


BLG208 


002. 


.4- 


-03. 


,7 


M 1-38 




18 


06 


05, 


.8 


-28 


40 


30 


T 


BLG128 


357, 


.3- 


-06, 


.5 


SB 50 




18 


06 


08, 


.2 


-34 


33 


30 


L 


BLG226 


003. 


.3- 


-03. 


.3 


PHR 1 806- 

X XXXV 1 UUv 


-2747 


18 


06 


28. 


,1 


-27 


47 


16 


L 


BLG217 


002. 


j. 


-03. 


1 


PHR 1806- 


-2824 


18 


06 


32, 


.2 


-28 


24 


27 


p 


BLG217,218 


002. 


.9- 


-03. 


.6 


MPA 1806 


-2812 


18 


06 


39. 


.3 


-28 


12 


21 


NULL 


BLG137 


357, 


.6- 


-06, 


.5 


PHR 1 806- 

X XXXV 1 Uvv 


-3416 


18 


06 


46, 


.3 


-34 


16 


04 


L 


BLG235 


003. 


.9- 


-03. 


,1 


KFL 7 




18 


06 


49, 


,9 


-27 


06 


19 


NULL 


BLG218 

J— 9 X— 4 VJ 1 U 


003, 


.0- 


-03. 


.6 


MPA 1 806 

-LVXX xX X \J\J\J 


-2807 


18 


06 


52, 


9 


-28 


07 


05 


L 


BLG185 

X-/ 1 > — J 1 


000, 
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